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Abstract 

A large Hilbert space is used for the calculation of the nuclear matrix ele- 
ments governing the light neutrino mass mediated mode of neutrinoless double 
beta decay (Oi^/3/3-decay) of "^^Ge, ^^^Mo, ^^^Cd, ^^^Te and ^^^Xe within the 
proton-neutron quasiparticle random phase approximation (pn-QRPA) and 
the renormalized QRPA with proton-neutron pairing (full-RQRPA) methods. 
We have found that the nuclear matrix elements obtained with the standard 
pn-QRPA for several nuclear transitions are extremely sensitive to the renor- 
malization of the particle-particle component of the residual interaction of 
the nuclear hamiltonian. Therefore the standard pn-QRPA does not guar- 
antee the necessary accuracy to allow us to extract a reliable limit on the 
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effective neutrino mass. This behaviour, aheady known from the calculation 
of the two-neutrino double beta decay matrix elements, manifests itself in 
the neutrinoless double-beta decay but only if a large model space is used. 
The full-RQRPA, which takes into account proton-neutron pairing and con- 
siders the Pauli principle in an approximate way, offers a stable solution in the 
physically acceptable region of the particle-particle strength. In this way more 
accurate values on the effective neutrino mass have been deduced from the 
experimental lower limits of the half-lifes of neutrinoless double beta decay. 
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Typeset using REVT^ 



I. INTRODUCTION 



The neutrinoless double beta decay (Oz//3/3-decay) continues to attract the attention of 
both experimentahsts and theoreticians for a long period. In this process the nucleus (A,Z) 
undergoes the transition to nucleus (A,Z+2) with the emission of two electrons: 

{A,Z) ^ {A,Z + 2) +2e~. (1) 

It is obvious that this process violates the lepton number Lg by two units and is forbid- 
den in the Standard Model. The background considerations imply that the Oi//5/?-decay is 
measured on nuclei for which the ordinary single beta decay is either forbidden by energy 
law conservation or strongly inhibited by large spin changes. The Oi//5/?-decay has not been 
seen in the experiment till now. A possible detection would be undoubtedly a signal of new 
physics beyond the standard model. The experimental lower limits provide us e.g. with the 
most stringent limits on the effective neutrino mass, and the parameters of right-handed 
currents and coupling constants of the supersymetric particles. If neutrinos turn out to be 
massive, the Oz//?/5-decay experiment is considered to be the most sensitive to the existence 
of Majorana neutrinos coupled to electron. It is wortwhile to notice that there is another 
two-neutrino mode of double beta (2i//5/3-decay) with two antineutrinos and electrons in 
the final state, which is allowed within the Standard Model. This mode being independent 
of the unknown particle physics parameters serves as a sensitive test of nuclear structure 
calculations. There exist extensive reviews of the theory and phenomenology of the double 
beta decay and we refer the interested reader e.g., to Refs. p for details. New contri- 
butions to the R-parity violating (-^p) supersymmetry (SUSY) mechanisms in Oz//5/5-decay 
are discussed in Ref. df. 

In order to deduce the limits on the parameters from the particle physics point of view, 
it is necessary to calculate the corresponding nuclear matrix elements. The proton-neutron 
Quasiparticle Random Phase Approximation (pn-QRPA) has been considered the most prac- 
tical method for nuclear structure calculations of nuclear systems which are far away from 
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closed shells [0. However, the extreme sensitivity of the calculated 2i//3/5-decay ma- 



trix elements in the physically acceptable region of the particle-particle strength of the 
nuclear Hamiltonian renters it difficult to make definite rate predictions This 



quenching behavior of the 2z//5/5-decay matrix elements is a puzzle and has attracted the 
attention of many theoreticians. Other shortcomings of the pn-QRPA (e.g., particle number 
non-conservation, the question of the proton-neutron pairing, the violation of the Pauli ex- 
clusion principle etc. ) indicate that we have to go beyond the pn-QRPA in the evaluation 
of the double beta decay nuclear matrix elements. Toivanen and Suhonen have proposed a 
proton-neutron renormalized QRPA (pn-RQRPA) to study the double beta decay |^ . The 
pn-QRPA is based on the renormalized quasiboson approximation, which considers the Pauli 
exclusion principle in an approximate way. Schwieger, Simkovic and Faessler have extended 
the pn-RQRPA to include proton-neutron pairing (full-RQRPA) |jl6[. The 2i//3/3-decay ma- 
trix elements calculated via pn-RQRPA and full-RQRPA have been found significantly less 



sensitive to the increasing strength of the particle-particle interaction |]T5[- [|T^. This fact 
indicated that both the Pauli exclusion principle and proton-neutron pairing play an im- 
portant role in the evaluation of the many-body Green functions for the double beta decay. 
In the meanwhile, some critical studies have shown that the renormalized QRPA has a few 
shortcomings, e.g., the violation of the Ikeda sum rule J. Hirsch et al. [Il8| , p!9| and 



J. Engel et al. ||2^ studied the validity of the renormalized QRPA within schematic exactly 
solvable models. Their studies confirmed that the renormalized QRPA offers advantages 
over the QRPA. But, they found some discrepancies between the exact and the RQRPA 
solutions after the point of collapse of the QRPA. However it is not clear whether their 
reults would also hold for realistic calculations with large model spaces and realistic effec- 
tive NN-interactions. In fact we do not know of any exactly solvable realistic model. The 



exact solution for the intermediate nuclear state discussed in Refs. |]T8| , p!9| within a schematic 
one level model space, after the point of collapse of QRPA, is below the initial ground state. 
This is not however the case of nuclei which undergo double beta decay . We note also that 
the violation of the Ikeda sum rule within the full-RQRPA with a large model space and a 



realistic effective NN-interaction is rather small (about 10-20 %) ||2T[ . 

We believe that the fuU-RQRPA is the most reliable method to deduce the desired 
interesting lepton number non-conserving parameters from the experimental lower limits 
on the half-lifes of neutrinoless double beta decay of heavier nuclei at present. To our 
knowledge the full-RQRPA has been applied for the first time by Simkovic et al. ||2T 
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calculations of Oz//?/?-decay. It has been found there by calculating the Oz//3/5-decay of ^^Ge 
via the pn-QRPA and full-RQRPA that by increasing the model space the pn-QRPA results 
are extremely sensitive to the renormalization of the particle-particle interaction in the 
physically acceptable region of the nuclear strength. This behaviour was similar to the 
one known from 2z//5/3-decay calculations. On the other hand the full-RQRPA results show 
an increased stability in respect to the stregth of the particle - particle interaction with 
increasing model space. This then sugests that the quenching of the 2z//3/?-decay matrix 
elements by the pn-QRPA is not a special phenomenon for the 2i//3/?-decay process but is 
a common feature to all many-body Green functions defining second order processes. As 
this quenching seems to have its origin in an inaccuracy of the calculation steming from 
the quasiboson approximation scheme it is necessary to recalculate the Oi//5/5-decay matrix 
elements with the full-RQRPA, which takes into account the Pauli exclusion principle and 
therefore one can expect to deduce more reliable limits on the lepton number non- conserving 
parameters . 

In this article, we shall study the nuclear matrix elements of light neutrino mass me- 
diated mode of Ouf3(3-decay of ^^Ge, i°°Mo, ^^^Cd, ^^^Te and ^^^Xe. The most stringent 
experimental lower limit of the life-time and large phase-space factors favour especially these 
processes for extracting an upper limit on the effective neutrino mass. Our first considera- 
tion is to find out whether the strong sensitivity of the Oi^/5/5-decay matrix elements within 
the pn-QRPA is a general feature to all double beta decay transitions. For that purpose we 
perform our calculations with a considerably large Hilbert model space . Our second task is 
to calculate the Oi//5/?-decay matrix elements within the full-RQRPA and investigate in this 
way whether one might be able to extract more accurate values on the effective neutrino 



mass. 



II. THEORY 

In the case of the hght neutrino mass mediated mode of the Oi//3/3-decay the weak beta 
decay Hamiltonian acquires the form: 

nf^ix) = ^2 [eL(x)7„i/e l(^)] ja{x) + h.c, (2) 

where ja{x) is the strangeness conserving charged hadron currents and ei,{x) and hix) are 
operators of the left components of fields of the electron and electron neutrino, respectively. 
We suppose that neutrino mixing does take place according to, 

^eL = E XkL, (3) 

k 

where, XkL is the field of the Majorana neutrinos with mass ruk and t/^ unitary mixing 
matrix. 

If we consider the usual approximations, i.e., non-relativistic momentum approximation 
for hadron currents and long-wave approximation, and replace the energies of the outgoing 
electrons in the denominators with the half of the available energy for this process, we get 
the following matrix element: 

' ' 2(27r)3 \^ VPl^ R 

u{p,){l-j,)Cu^{p2) (4) 

Here, pi, p2 are the four- vector momenta of the electrons and E\ are respectively the 
energies of the initial and final nuclear states. The effective neutrino mass is given as follows: 

3 

where exp(ici;j) is the CP eigenvalue of the neutrino mass eigenstate \xj >■ The nuclear 
matrix element of the process takes the form: 
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<a.. = M'J^ - (^^j M0^ (6) 

The Gamow- Teller and Fermi nuclear matrix elements take the form 

n. ^ /-"^ 9 < 0| II CrT(?^; J) II Jm >< Jm WO^^Uqr] L, J) \\ 0+ > 



mJ,L 







^ /o ' ^ * I * + - (E- + BO/2 1 ' ^ ' 



771 J 



with 



0°r(gr, J) = 5: 42V2 jz.(gr,) F;, (10) 

k 

Here, i? = tqA^/^ is the nuclear radius (ro = 1.1 fm), gy = 1.00, gA = 1-25 and qo ^ q for 
light neutrino. In the formulae of the and Mp" in Eqs. (|^) and @ it is somehow difficult 
to include the correlation function of the two interacting nucleons. Some attempts have been 
made by Krmpotic and Sharma [0. However, the operators Mq'j, and in and 
containing two one-body matrix elements are usually transformed to ones containing two- 
body matrix elements in relative coordinates by using the second quantization formalism. 
We then obtain: 

''°^\^Ur..r'^. (U) 

\ 1 / 



< O12 > = E i-r+'^'+'+^i2j + 1) 



jk ji J 

.^^r^,^ I ^ 

X <pk,pk'-J\f{n2)r^T^Oi2f{n2)\nl,nl'-J> 
X < 0/ II [Cpk'^ni']j II J'^mf >< .rmf\.rmi >< .Fmi \\ [c+ || 0+ > . (12) 

The short-range correlations between the two interacting nucleons are now taken into account 
by a correlation function P,|lO . 



l-e-""i2fi_5^2 



12; 



(13) 



with a = 1.1/m ^, b = 0.68/m ^. The neutrino-potential H{ri2) takes the form 



Hir) 



2R 



sin(gr) 



-dq. 



(14) 



nrJo g+(^]7.^ + ^]7/)/2(l + gVA2^4 
The parameter A of the dipole shape nucleon form factor is chosen to be 0.85 GeV [|5|,p!0 



Eq+ and il™/ = -Ejy 



For the half-hfe of the Oz//3/3-decay we obtain: 



K/2 



Oi/ 1-1 



mass) 



< ITLj, > 



rup 



(15) 



Goi is the integrated kinematical factor for the 0^^ — > 0^ transition [1T1,|3|,|23[| . 

In order to calculate the nuclear matrix element M^^^^ the full set of the intermediate 
nuclear states has to be constructed e.g., by the QRPA or RQRPA diagonalization. The 
fuU-RQRPA, which describes the excited states of the even-even nucleus, has been studied 



in Ref. |T6| and the pn-QRPA, which is a special case of the full-RQRPA in Refs. |[T4-[14 
Therefore, here we shall present only the formulae relevant to this work. 

The quasiparticle creation and annihilation operators (a^am^ and a^^^^, /x = 1, 2) are de- 
fined through the Hartree- Fock- Bogoliubov (HFB) transformation, which includes proton- 
proton, neutron- neutron and proton- neutron pairing p2| , p3| : 



/r+ \ 




( Uklp 


Uk2p 


—Vklp 






Ukln 


Uk2n 


— Vkln 






Vklp 


Vk2p 


Uklp 


\ Cnkrhk / 




\ Vkln 


Vk2n 


Ukln 



Uk2p 
Uk2n ) 



( (^tkmk \ 
^tknik 
^Ikni), 

V «2fcmfc ) 



(16) 



Here, c^ama i^rama) deuotcs the particle creation (annihilation) operator acting on a single 



particle level k with isospin t = p,n. The tilde ~ indicates the time reversed states Cr 



km,]. 



etc. 



In the full-RQRPA the commutator of two bifermion operators is replaced with its ex- 
pectation value in the correlated QRPA ground state \Qqrpa > (renormalized quasiboson 
approximation). We have 
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KrpaU,u{K /, J, M), {k\ I', J, M)]\0+^p^) 
1 1 



Ji 



Jk 



with jk = \^2jk + 1. The operator A^^^^k, I, J, M) creates a pair of quasiparticles coupled to 
angular momentum J with projection M. 



A%{k,l,J,M) = n{kiiM) E ^: 

n(fc/i, lu) = (1 + {-iy6ki6^,)/{l + 



(18) 



If we replace > iii Eq. (0) with the uncorrelated HFB ground state, we obtain the 

quasiboson approximation (i.e. 'D^k,vi\J'' ~ -'-)' which violates the Pauli exclusion principle 
by neglecting the terms coming from the commutator. The full-RQRPA takes into account 
the Pauli exclusion principle more carefully. The coefficients 'D^k,ui;j^i which renormalize 
the particle-hole and particle-particle interaction entering the A and B matrices of the full- 
RQRPA equation 

( I q\Ix'^\ 

-1 



( A k\ / X™ ^ 



A B 
B A 



Qm 

i 6 JTT 



V 



(19) 



are determined by solving numerically the system of non-linear equations [0,|T| 

T^k^lvJ^ 



1 ^ E '^kfik'tM'J'^'-J''^\^tMfi'i^^^'^ -J' 



3k fcV' 

1 



-2 E '^lvVv'J'^''^'^\^vy'{^t^'''^ 
Jl 



/ T/7r'M2 



(20) 



The selfconsistent scheme of the calculation of forward- (backward-) going free variational 
amplitude X (Y ), energies of the excited states Q'J't, and coefficients 'D^k,ui;j^ is a double 
iterative problem which requires the solution of coupled-non-linear equations. We note that 
in the limit T>^kM\J^ ~ ^ ^^"^ ^^^^ proton-neutron pairing is switched off, the solution 

of the full-RQRPA coincides with the solution of the pn-QRPA K 
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For the one-body densities in the full-RQRPA we can write: 



with m(/ia, vh) = . We note that the /, J"") and F^^(A;, I, J"") amphtudes 

are calculated by the renormalized QRPA equation only for the configurations fia < ub (i.e., 
/i = z/ and the orbitals are ordered a < b and fi = 1, u = 2 and the orbitals are not ordered) 
2|. For different configurations X^^{k,l,J'') and Y^^{k,l,J^) in Eqs. (|TD and (|D are 



(/) 



given as follows: 

x;Ak,i,r) = -{-iy'=+n-JX:^{i,k,r), (23) 

The index i (f) indicates that the quasiparticles and the excited states of the nucleus are 
defined with respect to the initial (final) nuclear ground state |0^ > (|0j >). We note that 
for Vk^i^j^ = 1 (i.e. there is no renormalization) and U2p = V2p = "Uin = Vin = (i.e. there 
is no proton- neutron pairing), Eqs. (0) and (p2D reduce to the expressions of the pn-QRPA 
0- 0. The overlap between two intermediate nuclear states belonging to two different sets 
is given by: 

<rmf\rm,>= E [T;j{ki,jnT;:{ki,r)-r;j{ki,jnr;:{ki,r)]. (25) 

fik<ul 



III. CALCULATION AND DISCUSSION 



We applied both the pn-QRPA and the full-RQRPA methods to calculate the Qvf3j3- 
decay of the A = 76, 100, 116, 128, 136 systems. In our calculations we tried to use as large 
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as possible Hilbert model spaces limited only by the power of the available computers. We 
have taken the folowing single particle model spaces for both, protons and neutrons : 

(i) For ™Ge Se, ^^^Mo Ru and "^Cd uiode\ space comprises 21 

levels: OS-^^j; OPl/2' 0^*3/2' -'-'^1/2' 0*^3/2' 0^5/2' ^Pl/2i ^Pzj2i O/5/2' O/7/2' '^^l/2i -'-'^3/2' -'-'^5/2' 
Ofi'7/2' 05'9/2' 2Pi/2' 2^3/2) 1/5/2' 1/7/2' 0^9/2' 0^11/2- 

(ii) For i28Te ^^^^ Xe and ^^^Xe ^^^^ we used 20 levels: 15^/2' 0^3/2, 0^5/2' lPi/2' 1^3/2' 

0/5/2' 0/7/2' 2S]^/2' 1*^3/2' -'-'^5/2' 05'7/2' 05'9/2' 2Pi/2' '^V'ij^i 1/5/2' I/7/2' 0^9/2' 0^11/2' 0^11/2' 
0^13/2- 

These model spaces are considerably larger as those used in any previous pn-QRPA calcu- 
lations 0- . The single particle energies have been calculated with a Coulomb-corrected 



Woods-Saxon potential. The nucleon-nucleon interaction used to calculate the nuclear wave 
functions is based on the Brueckner G-matrix derived from the Bonn one-boson-exchange 
potential, which in principle is a more consistent and better description of the NN-interaction 
in nuclei. The Brueckner reaction matrix is obtained by solving the Bethe-Goldstone equa- 
tion [E3]. The pn-QRPA is based on a Bardeen-Cooper-Schrieffer (BCS) transformation 



including proton-proton and neutron-neutron pairing correlations. In the case of the fuU- 
RQRPA the single quasiparticle energies and occupation amplitudes have been found by 
solving the HFB equation with p-n pairing in the above mentioned space Since our 



model spaces are finite, all pairing potential gaps are renormalized to the empirical gaps 
by the strength parameters (ipp, (i„„ a d^n- Technically this is achieved by performing a 
BCS and HFB calculation and comparing the obtained pairing gaps with the ones extracted 



from the empirical separation energies in a manner described in Ref. |2^. The renormaliza- 
tion parameters (ipp, dnn and d^n together with the experimental proton ( Ap^^ ), neutron 
( A^'^P ) and proton - neutron ( 5^^^ ) pairing gaps, for all studied nuclear systems, are 



listed in Table |. The experimental pairing gaps are defined by Moeller and Nix p5| , |26 
By glancing at the Table | we see that the d^^ and (i„„ values are close to unity and d^n 
is higher than these values. It is because for spherical nuclei the HFB-transformation can 
only describe correlations for pairs with J=0 and T=l and not for pairs with J=0 and T=0. 
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The T=0 pairing is effectively taken into account by the renormahzation of the T=l J=0 
n-p interaction leading to a higher value of dp^. However, we do not want to focuse our 
attention to the problem of proton-neutron pairing, which seems to play a significant role 



in the QRPA calculation of the double beta decay process |T6| , |23| but it is less important in 
the case of the renormalized QRPA calculations |2T[]. The proton-neutron pairing problem 
is extensively discussed in p2| , ^ and references cited there. We note that in the BCS and 
HFB calculations we neglected the mixing of different "n" but the same "Ijm" orbitals. We 
suppose that shell mixing is not significantly affecting the BCS and HFB solutions because 
their off-diagonal pairing matrix elements are quite small. 

After settling the values of the pairing parameters, the parameters which remain to be 
fixed are the particle- particle and the particle- hole strengths. The particle - particle and 
particle - hole channels of the G-matrix interaction are renormalized by introducing the two 
parameters and gpf^^, which, in principle, should be close to unity. Our adopted values 



were g^h = 0.8, as in our previous calculations [£1|)E31 9pp varied in the interval 0.70 - 
1.30 which can be regarded as physical. 

The nuclear matrix elements M^^^^^ for the most interesting nuclei obtained within the 
pn-QRPA are shown in Fig. 1 (a). We see that M^^^^ for the A=76, 100, 116, 128 and 136 
systems, becomes singular with increasing strength of the particle-particle interaction and 
even crosses zero in the physically acceptable region of the parameter gpp . This behavior 



has not been found in the previous calculations 0- |T0[ because the model spaces used there. 



were too small. To our knowledge only Krmpotic and Sharma |Tl|] have studied the model 
space dependence of the Oi//5/?-decay matrix elements. For a relatively large model space 
they found for the Oz//3/3-decay of ^*Ca a similar behavior. However this was further out 
of the physicall region of the strength renormahzation parameters. We note that there is 
a principal difference between our calculations and those of Krmpotic and Sharma. They 
used a zero range 5-force interaction and a different treatment to the two-nucleon correlation 
function, which were incorporated in the formalism in Eqs. (^0). To our opinion their type 
of two-nucleon correlations influences the result unsignificantly. The two-nucleon correlations 
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presented in this work give an effect of about 30-40%. We should note further that ^^Ca is 
a closed shell nucleus and therefore it is not very suitable for the QRPA calculation. 

The QRPA quenching mechanism for M^^^^ in Fig. 1 (a) has its origin, exactly in 
the same way as the quenching mechanism of 2i//5/5-decay matrix elements, i.e., in the 
generation of too much ground state correlations with increasing Qpp near the collapse of 
the QRPA. A larger model space means more ground state correlations, i.e., a collapse 
of the QRPA solution for smaller gpp. As a consequence the validity of the quasiboson 
approximation in the evaluation of the Oz//3/3-decay matrix elements is questioned because 
of the generation of too much ground state correlations. For that reason it is necessary to 
perform the calculation of M^^^^ in the framework of the renormalized QRPA, which takes 
into account the Pauli exclusion principle in an approximate way and also considers proton- 
neutron pairing correlations. In Fig. 1 (b) we present our results with the full-RQRPA 
method. From the comparison with the Fig. 1 (a) it follows that the inclusion of ground 
state correlations beyond the QRPA in the calculation of M^^^^ removes the difficulties 
associated with the extreme sensitivity of M^^^^ on the particle-particle strength. The 
strong differences between the results of both methods indicate that the Pauli exclusion 
principle plays an important role in the evaluation of the 0z//3/3-decay. 

A weakly dependence of the M^^^^^ on gpp inside the physical range 0.8 < gpp < 1.2 
allows us to have more confidence for deducing the effective neutrino mass < rrii, > from 
the available experimental lower limits on the half-lives of the Oi//3/?-decays Ty^'^^^. The 
nuclear matrix elements M^'^^g obtained within the full-RQRPA for gpp =1.0, the integrated 
kinematical factors Gqi, T^^'^^^ and the limits on the effective neutrino mass | < niy > \ 
deduced from M^'^^g {gpp =1.0) and Ty^'^^^ are listed in Table ||. A more stringent upper 
limit on < rrij^ > is favored by large values of M^^^^, Gqi and T^^'^^^. We see that the most 
stringent upper limit < rrii, > < 1.1 eV is deduced for the A=76 system mainly because of 
the unbelievable high upper limit on T^^^^^ given by the Heidelberg- Moscow collaboration 



1 2^. Also of interest is the value T'^J^ ^'^^ (see Table |I|), which is calculated by assuming 



< nil, > I = leV. This value indicates that for further experimental measurements the 
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most perspective candidate is ^°°Mo. 

IV. SUMMARY AND CONCLUSION 

In summary, we have studied the nuclear matrix elements entering the light neutrino 
mediated mode of the neutrinoless double beta decay of some experimentally interesting 
nuclear systems, A=76, 100, 116, 128 and 136 . The calculations have been performed within 
both, the pn-QRPA and the fuU-RQRPA methods with a large Hilbert model space. We 
have found that in the framework of the pn-QRPA for a large enough model space the OuPP- 
decay matrix elements demonstrate an instability with respect to the renormalization of the 
particle-particle strength, similar to the one known from the 2i//3/3-decay mode. The value of 
the matrix element crosses zero and it is then difficult to make definite rate predictions. We 
believe that the common quenching phenomenon which is independent of the studied nucleus 
and double beta decay process could have its origin only in the approximation scheme. The 
fuU-RQRPA which includes proton-neutron pairing and the Pauli effect of fermion pairs goes 
beyond the quasi-boson approximation. The inclusion of the Pauli principle eliminates the 
instabilities that plaque the pn-QRPA. The Oi//3/3-decay matrix elements calculated via the 
fuU-RQRPA are stable in respect to the changes of the particle-particle force and it allows 
us to deduce more accurate limits on the effective neutrino mass. The largest Oi//3/3-decay 
matrix elements are 4.22 and 3.28 associated with the A = 100 and 128 systems, respectively. 
A large value of the matrix element together with a large value of the kinematical factor 
favour especially ^°°Mo for further experimental study of the 0i//3/5-decay. At present the 
A= 76 and 128 systems provide us with the most stringent limit on the effective neutrino 
mass i.e., 1.1 - 1.2 eV. 
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TABLES 

TABLE I. Experimental proton ( A^^^ ) , neutron ( A^^^ ) and proton - neutron ( ) pairing 
gaps and renormalization constants of the proton - proton ( dpp ) , neutron - neutron ( dnn ) and 
proton - neutron ( dpn ) pairing interactions for all nuclei studied here. 



Nucleus 


( A<^^'P A'^^'P A'^^P 
( [MeV] ) 


model 
space 


dpp 


dnn 


dpn 




(1.561, 1.535, 0.388) 


21 level 


0.899 


1.028 


1.506 


115642 


(1.751, 1.710, 0.459) 


21 level 


0.934 


1.059 


1.325 




(1.612, 1.358, 0.635) 


21 level 


0.980 


0.923 


1.766 




(1.548, 1.296, 0.277) 


21 level 


1.002 


0.945 


1.568 


ll^Cdes 


(1.493, 1.377, 0.371) 


21 level 


0.953 


0.922 


1.822 


50^ C dee 


(1.763, 1.204, 0.128) 


21 level 


1.00 


0.873 


1.460 




(1.127, 1.177, 0.149) 


20 level 


0.873 


0.942 


1.780 


ifXe74 


(1.307, 1.266, 0.199) 


20 level 


0.920 


0.972 


1.530 


Ae82 


(0.971, 1.408, 0.0) 


20 level 


0.771 


0.803 


0.0 




(1.245, 1.032, 0.165) 


20 level 


0.875 


0.899 


1.716 
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TABLE II. The nuclear matrix elements M^^^^ (see Eqs. (6) and (11-12)) obtained within 
the full-RQRPA for gpp = 1.0, the integrated kinematical factors Gqi, the limits on the effective 
neutrino mass | < rrii, > | deduced from the experimental limit of the Oz^/3/3-decay lifetime T-Pjl"^^^ 



for the nuclei studied in this work. ^ calculated Oi^/5/5-decay half-life times assuming 

I < ray > I = leV. 



nucleus 


mass 


( rm 

^ Ul 

J-|^g-14 ygQT-g-l] 


rpOu-leV 
1/2 

[years] 


rpOu—exp 
1/2 

[years] ref. 




1 < Uli, > 

[eV] 




1.86 


0.7928 


9.5 X 10^4 


> 7.4 X 1024 (90% C.L.) 


m 


<1.1 




4.22 


5.731 


2.6 X 10^3 


> 4.4 X 1022 (68% C.L.) 


m 


<2.4 




2.47 


6.237 


6.9 X 10^3 


> 2.9 X 1022 (90% C.L.) 


m 


< 4.9 




3.28 


0.2207 


1.1 X 10^5 


> 7.3 X 1024 (68% C.L.) 


m 


<1.2 


i36xe 


0.96 


5.914 


3.4 X 10^3 


> 6.4 X 1023 (90% C.L.) 


[32 


< 3.7 
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FIGURES 



FIG. 1. The calculated nuclear matrix element M^^^^ for the 0z//3/?-decay of ^^Ge, ^^^Mo, 
^^^Cd, ^^^Te and ^^^Xe as a function of the particle-particle interaction strength gpp. In (a) M^^^^^ 
has been calculated with the pn-QRPA method, in (b) with the full-RQRPA . 



19 



FIGURES 





